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Abstract 

The cross sections for the formation of scalar resonances, leptoquarks or 
squarks, in electron/positron-proton collisions at HERA are presented includ- 
ing next-to-leading order QCD corrections. Depending mildly on the mass 
of the resonances, the ^-factors increase the production cross sections by up 
to 30% if the target quarks are valence quarks. The QCD corrections to the 
partial decay widths of leptoquarks/squarks to leptons and quarks are small. 
The electron spectrum in the decays is softened nevertheless by perturbative 
gluon radiation at a level of 3.4 GeV for a leptoquark/squark mass of 200 
GeV. 
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1 Introduction 



An excess of events has recently been observed in deep-inelastic e + p scattering at HERA 
for high x and high Q 2 |], |2|. If these events are attributed to the formation of a resonance 
in e + q collisions, decaying subsequently into lepton-quark final states, the events should 
cluster at the invariant mass M = \/xs. Moreover, if scalar resonances are formed, the 
distribution in the scaling variable y should be flat, corresponding to isotropic angular 
decay distributions in the rest frame of the resonance. Indeed, a clustering around M = 
200 GeV is indicated in the HI data and the events are distributed up to large y values 
|Ij] . These observations have initiated theoretical analyses of possible resonance formation 
at HERA, and cross checks with experimental results from LEP ||, the Tevatron Q and 
from rare decays have been performed. The analyses have been carried out for leptoquarks 
in general and squarks in supersymmetric theories with i?-parity breaking in particular 

In the present report we present the cross sections for the formation of scalar reso- 
nances, leptoquarks or squarks, including next-to-leading order QCD corrections in e ±r p 
collisions. The step beyond the leading order analysis [|TT|] is motivated by two points. 
First, the theoretical predictions for the cross sections are stabilized with regard to varia- 
tions of the renormalization and factorization scales. The choice of the leptoquark/squark 
mass M, which can a priori be considered as a natural choice for these spurious param- 
eters, can thus be investigated in detail. Second, to derive the theoretical predictions for 
the formation cross sections to an accuracy of 0(as) ~ 10% is a reasonable goal for the 
time being. The results are presented for a single scalar resonance 

e + q -> IQ/q (1) 

where e generically denotes electrons and positrons, q quarks and antiquarks, and IQ/q 
leptoquarks and squarks (Fig. |l]a). The coupling A is defined by 

tint = A(e q)S + h.c. (2) 

with the (opposite) helicities of the lepton and the quark fixed for the specific lepto- 
quark/squark state S (see e.g. Table 1 of Ref. 0). 

By the same token, we have determined the partial width of leptoquark/squark decays 

IQ/q - e + q(g) (3) 

to order a s in QCD. For massless quarks in the final state, the (scaled) lepton energy 
x e = 2E e /M is fixed to be 1 in lowest order. The lepton energy is softened by gluon 
radiation and the spectrum becomes continuous, peaking however strongly for x e just 
below 1. 
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2 Production of Scalar Particles in ep Collisions 



The cross section in leading order QCD can be cast into the form 

o~lo 
0o 



ao x q (x, M 2 ) 
vrA 2 



4M 2 



(4) 



with 



x 



M 2 /s (5) 

The density of the target quarks at the scale /ijr = M is denoted by q(x, M 2 ). 

In addition to the standard mass corrections for colored quarks and bosons which 
are defined on the mass shells, the QCD corrections involve the renormalization of the 
e — q — LQ/q vertex displayed in Fig. [l|b. The e — q — LQ/q coupling A has been defined 
at the point = M in the MS renormalization scheme. The QCD vertex corrections 
make the coupling run according to the formula 



A 2 M 



A 2 (M) 



(6) 



The second class of diagrams, Fig. |T]c/ d, describe the radiation of gluons off the initial 
quark and final leptoquark/squark state: e + q — > g + LQ/q. The sum of the virtual 
and bremsstrahlung corrections is ultraviolet and infrared finite. Finally, the formation 
process initiated by lepton-gluon collisions must be added: e + g — > q + IQ/q, as shown in 
Fig. [TJe. The collinear divergences due to gluon emission off the initial state quark line and 
due to gluon splitting to quark-antiquark pairs will be absorbed in the renormalization of 
the quark and gluon parton densities [|12| . 



The final result of this simple analysis can be summarized in the standard form 



a(ep -> LQ/q + X) 



0-Q 



2tt 2 a s 



xq(x, M 2 ) + Aa q + Aa g 



(7) 



with 



rl Ct ( Z 4 

Aa q = / dx' q(x', M ) x -^a --P qq {z) \ogz + - [1 + z 

/ log(l — Z 



+ - 



1 - z 



1 - Z 



(2 + z + z 2 ) log(l - z) 



rl Ct Z \ 

Aa g = / dx' g{x',M 2 ) x ^c-o- \-P qg {z) log 

Jx 7T I y 
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+ z(l — z) log z + 1 
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The auxiliary variable z is defined as usual by z = x/x' . P qq and P qg are the standard 
q — > q(g) and g — > q(q) splitting functions 



The distribution F + is defined by F + (z) = F(z) — 5(1 — z) J dz'F(z'), regularizing the 
singularities at z — 1. The renormalization point hr and the factorization point fip have 
been identified with the mass M of the resonance. They may be re-introduced as free 
parameters in Eq. (f|) by substituting A(M) — > \(/j,r) and adding the term + log ji 2 R /M 2 
to the coefficient of as/rc in the virtual correction; in addition, the mass scale M in the 
parton densities must be replaced by fip and log \ipjM 2 must be added to the logarithms 
attached to the parton splitting functions. 

The results are exemplified in Fig. |2| for the .fT-factors properly defined by 



where all quantities (couplings, parton densities) in the numerator and denominator are 
to be evaluated consistently in next-to-leading and leading order, respectively. This defi- 
nition guarantees the correct theoretical interpretation of the size of the QCD corrections. 
The i^-factors are presented for e ± projectiles on valence (u, d) targets in Fig. |2|. The 
.fT-factors are larger than unity in the range of leptoquark masses we are interested in. 
The NLO cross sections are larger than the leading order cross sections by about 20 to 
30%. The vertex corrections reduce the cross section, as apparent from Eq. (|^). This 
effect is over-compensated however by the quark contribution Aa q , while gluon splitting, 
Acr g , plays a minor role. For sea quark targets, the .fT-factor is close to 1.15 for M = 150 
GeV and falls to values below 1 for masses above 195 GeV. The different behavior of 
if-factors for valence and sea quarks is a straight consequence of the different shapes of 
the parton densities. 

The improvement of the theoretical prediction is clearly visible in Fig. [3]. The cross 
sections in next-to-leading order are compared with the leading order for the process 
e + d — > IQ/q as a function of the renormalization scale \ir and the factorization scale \if- 
Hr and \ip are identified for the sake of simplicity. While the LO cross section changes by 
nearly a factor of 2 if /x/M is varied from 3 to 1/3, the change of the NLO cross section is 
damped to 1.2 in this range. Thus, if the next-to-leading order corrections are taken into 
account, the theoretical prediction of the cross section for scalar resonance formation is 
stabilized considerably. 




(9) 



8 



K = cr NLO /a LO 



(10) 



4 



3 Decay Width and Lepton Spectrum 



The partial width for the decays 

IQ/q^e + q(g) (11) 
can easily be calculated in a similar way including the one-loop QCD corrections: 



T[IQ/q^eq{g)\ 



X 2 M 

167T 




(12) 



The coefficient in front of as/ft is tiny, ~ 0.11, so that QCD corrections to the width 
are small. For small Yukawa couplings A ~ e/10 the partial width is therefore very 
narrow, i.e. T ~ 3 MeV for M ~ 200 GeV. If for leptoquarks/squarks no other channels 
are open, apart from LQ — > eq or uq, and q — *> eq, the lifetime is very long, and the 
lepto quark/ squark will be bound into a fermionic {LQq) or (qq) state. This is different for 
squarks if the branching ratio for i?-parity violating decays is small compared to i?-parity 
conserving (neutralino/gaugino) decays. 

To lowest order, the scaled lepton energy of the decay process (|]) in the rest-frame[] 
of the isolated leptoquark or squark, 



x 



2E e /M with < x e < 1 



(13) 



is fixed to 1 if the final-state quark is treated as a massless parton. However, perturbative 
gluon radiation (and the non-perturbative hadronization of quark jets) reduces the lepton 
energy and gives rise to a continuous energy spectrum. This effect is analogous to the 
radiation of gluons in qq pair production of e + e~ annihilation; LEP measurements can 
therefore serve as a solid basis for estimating non-perturbative parameters which are not 
accessible to theoretical QCD analyses. 

The energy spectrum of the lepton, after real gluon emission is switched on, is given 
by the expression 



x, < 1 : 



1 dY 



r dx e 



pQCD 
0(a s ) 



4a, 

3 7T 



1 / 1 7\ , 1 7 x e 
— log + log — + - - — 

Xx V Xi 4/ X\ 4 4 . 



(14) 



with the abbreviation x\ = 1 — x e . This form of the spectrum is normalized to the partial 
eq decay width. (If other decay channels are open, Eq. ( |14"D must be supplemented by 
the branching ratio B eq for absolute normalization). In the limit x e — > 1, the leading 
singularities may be resummed to give 



x ( - 



1 : 



1 dY 


pQCD 


1 dF 


pQCD 

x exp 

0(a s ) 




■ 1 


r dx e 


resum 
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- - log — 

4 xi 



(15) 



lr rhe scaled lepton energy can be defined in a Lorentz-invariant way as x e — 2{p e PLQ/q)/M 2 
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The singularity of Eq. fll4] ) for x e — > 1 is removed by resummation |14| |. In fact, the 



spectrum flT5| ) approaches zero for x e — > 1 and a sharp maximum develops just below 
1. The final lepton-energy spectrum as predicted in pQCD is shown by the full line in 
Fig. The spectrum in the approximation of one-gluon radiation without resummation 
(singular near x e — > 1) is indicated by the broken line. 

The softening of the lepton spectrum due to the perturbative gluon radiation can be 
characterized by the average lepton energy (x e ) defined for LQ — > e + q(g) decays. It is 
easy to derive from the spectrum: 

(x e y QCD = i-^ (i6) 

9 7T 

This corresponds to a loss of ~ 3.4 GeV in the lepton energy due to perturbative gluon 
radiation for a leptoquark mass of 200 GeV. 

In addition to the perturbative energy loss, residual non-perturbative hadronization 
effects reduce the lepton energy, affecting the shape as well as the average energy. Such 

effects can be parametrized in a power series 5x? p = -d/M H . Estimates of the 

size of the non-perturbative coefficient C\ may be derived from the shape variable thrust 



15 1 . A closer analogon is the direct-photon spectrum in e + e~ annihilation to hadrons 
|I73fl at LEP. With C\ ~ 0.7 GeV for thrust, the non-perturbative contribution Sx^ p ~ 
(— a few 10~ 3 ) is expected well below the perturbative attenuation effect. In Monte 
Carlo calculations, which are beyond the scope of the present analysis, this qualitative 
estimate can be turned into a quantitative prediction. The overall picture drawn from the 
pQCD calculation presented above coincides well with the result of shower Monte Carlo 
programs employed in the HI analysis fll? |. This is expected from the QCD analyses at 



LEP. Scenarios in which squarks decay so fast that the color flux lines of the quarks in 
the decay final state are connected directly to the proton remnant, require more elaborate 
Monte Carlo programs. 
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Figure 1: The basic diagrams contributing in QCD to the formation of lepto- 
quarks/s quarks in ep collisions, (a) Born term; (b) vertex loop; (c) gluon radiation off 
the leptoquark/squark; (d) gluon radiation off the quark; (e) gluon splitting to quarks and 
antiquarks. The same diagrams (a) to (d), interpreted backwards, describe the decay of 
leptoquarks/s quarks to order a s . 
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Figure 2: K -factors for ed, eu — > LQ /q as a function of the leptoquark/s quark mass. 
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Figure 3: Comparison of the renormalization and factorization scale dependence in LO 
and NLO for the cross section for a(e + d — > IQ/q). 
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Figure 4: TTie lepton-energy spectrum as predicted in perturbative QCD (full line). (The 
broken line corresponds to the 0(a s ) result without resummation near x e — > 1). 
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